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Abstract 

Conformal field theory and Bethe ansatz are used to investigate the low energy 

features of the spectral function in one dimensional models which exhibit a 

gap in the spin or in the charge excitation spectrum. Exotic behavior is found 

in the superconducting case, where the Green function displays momentum 

dependent Luttinger Liquid exponents. The predictions of the formalism are 

confirmed by Lanczos diagonalizations in the t J model up to 32 sites. These 

results may be relevant in connection to photoemission experiments in quasi 

one dimensional insulators or superconductors. 
75.10.Jm,74.25.Gz,71.27.+a 
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In the last few years direct and inverse photoemission experiments have considerably 
improved, allowing for the accurate determination of momentum dependent energy spectra 
in low dimensional systems • A first principle interpretation of this class of experiments 
requires a deep understanding of the effects of correlations on the electron (or hole) spectral 
function. One dimensional (ID) metals [Q], have been the subject of an intense theoretical 
effort which led to a complete characterization of the long wavelength, low energy properties 
of electron dynamics and to the actual calculation of the correlation exponents which appear 
in the electron Green function However, some doubt has been cast on the relevance 

of these calculations for photoemission experiments because most of quasi ID systems are 
close to density wave or superconducting instabilities which open a gap in the charge or 
spin spectrum ||. If the excitation spectrum is fully gapped, the Green function takes a 
free particle-like form in every dimension. Instead, the effects of a a single branch of gapless 
excitations (either spin or charge) have not been addressed in detail before, probably because 
the prejudice prevailed that a system with a gap should display exponentially decaying 
correlation functions both in space and time. 

In this Letter, we develop a microscopic theory for determining the low energy properties 
of the spectral function in one dimensional correlated electron models with a gap either in 
the charge or in the spin channel. The main result of this work is that the presence of gapless 
excitations induces anomalous exponents in the Green function, by a non trivial interaction 
with the extra electron (or hole) injected into the system. As a consequence, we generally 
find a spectral function with singularities along lines in the (k, u) plane. These singularities 
are characterized by critical exponents which possibly depend on the momentum k of the 
electron: Numerical diagonalizations in the t J model at J = 2t fully confirm this picture 
providing quantitative agreement with predictions based on conformal field theory and Bethe 
ansatz techniques. 

Spin and charge gaps are treated on the same footing in this Letter because they are 
believed to give rise to the same kind of singularities in the spectral function, despite the 
quite distinct physical nature of the state. In fact, it is possible to build specific models 
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where the two regimes are mapped one onto the other: A well known example is the negative 
U Hubbard model at zero magnetization which is the prototype of one dimensional "super- 
conductors" i.e. ID models with spin gap and quasi off diagonal long range order. This 
model, via a particle-hole transformation, is mapped into the half filled, positive U system 
with non vanishing magnetization which, on the contrary, is a Mott insulator characterized 
by a charge gap. The full Green function remains unaltered by particle-hole transformation 
leading to the same photoemission spectra in the two cases. In the following, we will ex- 
plicitly deal only with the repulsive, half filled case, at arbitrary magnetization. The results 
will however hold both for one dimensional insulators and "superconductors" , being related 
only to the presence of a branch of gapless excitations in the spectrum. 

The quantity which we are going to investigate is the hole Green function: 

G(p,t) = i < tt | cl a e-^ H - E °-^ c p , CT |* > 6{t) (1) 

where |*q > (E Q ) is the ground state (energy) of the system with no holes and 9(f) is the step 
function. Due to spin-charge decoupling, the total energy E and momentum p are naturally 
written as a sum of a holon term e^ik) and a spinon contribution e s (Q) with p = k + Q. At 
long wavelengths it is known that holons and spinons behave as independent particles whose 
dynamics is governed by two commuting hamiltonians: H c and H s respectively. By substi- 
tuting this decomposition H = H c + H s into Eq. ([!]) and taking momentum conservation 
into account, we find that the hole Green function can be written as a sort of convolution 
between a holon (Gh) and a spinon (Z) term: 

G( P ,t) = J^ Gh (p-Q,t)Z p (Q,t) (2) 

where Gh(k, t) is just a free propagator: Im Gh(k, u) = n 5{uj—e h (k)). This simple form of the 
holon Green function is due to the presence of a gap in the charge excitation spectrum of the 
model which confines the low energy processes in the single holon sector. Instead, the spinon 
contribution Z is highly non-trivial due to the presence of gapless excitations and exhibits 
anomalous exponents at particular momenta Q v . As a result, the most relevant singularity 
in the full spectral function occurs at frequencies determined by the hole dispersion: 
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Here X v {p) is a momentum dependent critical exponent determined by the low energy prop- 
erties of the spinon dynamics, while Q u = (2z/+ 1)Qf is an odd multiple of the spinon Fermi 
momentum Qp = 7r(l/2 — m), m being the magnetization per site. Equation (||) generalizes 
the exact form found in the U — > oo limit of the half filled Hubbard model @ where Z P (Q, t) 
was explicitly calculated and turned out to be independent of p and t. In this limit, the hole 
dispersion is th{k) = — 2cos& and the critical exponent takes the value Xo = X-\ = 1/4. 

In the following, we will analyze the long wavelength behavior of Z p (Q,t) which leads 
to Eq. (|3|) in the particular case of single (spin down) hole in the t J model at arbitrary 
magnetization m. The case m = gives information about the Mott insulator while the 
m > (m < 0) choice refers to photoemission (inverse photoemission) experiments in 
"superconductors" via spin-up (spin-down) particle-hole transformation in the less than 
half filled attractive Hubbard model. Due to the universality underlying the behavior of one 
dimensional physics, we expect that these results will be qualitatively valid for generic ID 
electron system displaying a gap in the excitation spectrum. In fact, it is well known || 
that in this case the renormalization group flow drives the model towards the Luther Emery 
fixed point irrespective of the details of the microscopic hamiltonian. On the other hand, 
the tJ model allows for a direct comparison with Lanczos diagonalizations which can be 
pushed to fairly large lattice size in such a system. 

The single hole problem in the t J model can be reduced to a pure spin problem by a 
Galileo transformation which fixes the hole at the origin O of the L-site lattice (L is 
chosen to be even). In this way, the charge degree of freedom can be exactly traced out 
leaving the problem of an effective momentum dependent spin hamiltonian: 



Here p is the total lattice momentum of the one hole state and T is the translation operator 
along the squeezed chain with I = L — 1 sites, i.e. without the origin O. The hamiltonian 



L-2 




e ip T + e - ip T^] + J ]T Si ■ S 



(4) 



i=l 



4 



H e H is written as the sum of two physically different contributions: The magnetic term is 
the standard Heisenberg model with open boundary conditions because no spin is present 
at the origin, while the kinetic of the hole manifests itself via the action of the translation 
operator T. 

This mapping of the hole problem into a spin hamiltonian is exact for every J and allows 
to interpret the presence of the hole as the inclusion of a special type of boundary operator in 
the bulk spin hamiltonian. An insight to the general features of the energy spectrum of this 
hamiltonian can be obtained by examining the J — > limit, where all the spin configurations 
on the squeezed chain which are eigenstates of the translation operator T\ipQ >= e t( ^\tpQ > 
are degenerate provided they correspond to the same (spinon) momentum Q. This degener- 
acy is lifted by the magnetic term, which, at first order in J, selects the lowest energy state 
of the Heisenberg ring with the given spinon momentum Q. The corresponding hole energy 
is E = th(p — Q) + e s (Q), where the holon band is €h(k) = —2 cos A; and the O(J) spinon 
dispersion e s (Q) only depends on the bulk properties of the Heisenberg model. In this limit, 
the effects of spin-charge decoupling on the energy spectrum come out rather naturally as 
well as the role of the hole kinetic contribution in modifying the boundary conditions of the 
Heisenberg model, from open to periodic. Due to the peculiar form of the hole boundary 
operator in Eq. (Q), the long wavelength behavior of H e ^ is associated to a new class of 



fixed points different from those found in the framework of the static impurity problem |10 

After a standard Jordan- Wigner transformation the spin hamiltonian maps onto an elec- 
tron system of interacting spinless fermions at density p = | — m. At low energy the relevant 
degrees of freedom for a many fermion system are the momenta close to the Fermi points 
±Qf- It is then possible to take the continuum limit of the model defining two independent 
fermionic fields iPr(x) and for the right k ~ Qp and left k ~ —Qf movers on the 

squeezed chain with < x < I. The long wavelength limit of the translation operator can 
be written in terms of a well defined spinon momentum operator P: T = e lP where 



P = Qf /<fc{K(*)«*)-<*t(x) W z) 



+ 



i[^ R (x)d x ip R (x) + ipl(x)d x tl) L (x)\ } (5) 

As anticipated, the competition between the hole kinetic term and the magnetic interaction 
in H e ff selects a particular effective boundary condition for the interacting fermion gas which 
simulates the presence of a scattering potential at the boundary together with a magnetic 
flux across the ring. This gives rise to independent boundary conditions for the right and 
left movers defined by two phase shifts: 

4>Ux + l) = e lS ^i(x) (6) 

The long wavelength analysis proceeds by noting that, due to spin charge decoupling, the 
effective hamiltonian splits into the sum of two commuting terms 

H eff = t h (p-P) + Hj (7) 

governing charge and spin dynamics respectively. Here Hj is the usual long wavelength 
form of the Heisenberg hamiltonian in terms of the fermionic fields (ip^, ip R ) M with the 
particular boundary conditions ©. It is known that the interacting hamiltonian (0) can be 
turned into a free Fermi problem by a canonical transformation which further modifies 
the boundary conditions (@). The bulk Luttinger liquid properties of Hj are described by 
the dressed charge K p shown in Fig. 1 ||. 

Having characterized the long wavelength physics of the effective hole hamiltonian (|J) 
in terms of the dressed charge and the two phase shifts, let us move to the analysis of the 
spectral properties of the hole motion. As a first step, we fix our attention on the overlap ( p 
between the pure Heisenberg ground state on L sites j^o > and the one hole ground state 
|\&p > at fixed momentum p: 

( P =\< %\c P ,i\^o > | 2 oc L- 2X ^ (8) 

This quantity naturally enters the calculation of the hole spectral function as can be im- 
mediately checked by use of Lehmann representation. As indicated in Eq. (jS)), the overlap 
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(p vanishes in the thermodynamic limit with a critical exponent which can be explicitly 
evaluated in terms of the previously introduced phase shifts: 

x ^ Kf (^ + v J + ^(^J (9) 

The calculation parallels the known derivation of the orthogonality catastrophe in the im- 
purity problem ||11|| . Here, the integer number v defines the total spinon momentum of the 
intermediate states Q u = {2u+1)Qf corresponding to an odd number of low energy spinons. 

The phase shifts (||) are non universal depending on the short wavelength properties 
of the model, however it is possible to relate them to the form of the energy spectrum 
E(p) of one hole at fixed momentum p which, according to Eq. ([j]), can be written as 
a sum of a charge and a spin contribution. In fact, we expect that low energy spinon 
dynamics is governed by some effective conformal field theory which should reflect on the 
structure of the size corrections to the spinon contribution to the ground state energy: 
E{p) - Le^p) = (AE c (p) + AE s {p)) /L with 

AE s {p) = -v. \ + 2nv s X u (p) (10) 
o 

Here v s is the spinon velocity: v s = de s (Q)/dQ evaluated at Qf and X v {p) coincides with 
the critical exponent (||). 

In order to determine the unknown phase shifts 5r and 5l, we have analyzed the size 
corrections to the ground state energy (at fixed momentum p) in the Bethe ansatz soluble 
models: The Hubbard model at U > |12| and the t J model at J = 2 |13[ with one hole 



and arbitrary magnetization. By suitably generalizing the pioneering work of Woynarovich 
14| to the single hole case, we found exactly the form (|T0|) of the energy size corrections, 



with quantitative predictions for the phase shifts which in fact explicitly depend on the total 
momentum p at every non-zero magnetization m. For m — 0, i.e. for the Mott insulator 
case, instead, we always find that only one of the two phase shifts is different from zero 
and takes the value tt, both in the Hubbard and in the tJ model. Another analytic limit 
is the J — >• at arbitrary magnetization m where again the phase shifts are independent of 



p but are functions of the magnetization: 6r = tt (1 — m) and 5l = nm. Figure 2 shows 
the exponent X (p) as a function of the total momentum p for the Bethe ansatz solvable 
limit of the t J model (J = 2) at several magnetizations. A comparison with the value of 
the overlap exponent obtained by use of Eq. (|8]) through Lanczos diagonalization of the 
model with magnetization m = ±0.25 is shown in Fig. 3 . Lattice sizes ranging from 16 
up to 32 sites have been used to fit the exponent 2X u (p) in Eq. (||) leading to a quite 
good agreement between analytical and numerical results. This comparison gives confidence 
on the interpretation of the Bethe ansatz results for the single hole size correction in the 
framework of conformal field theory. 

Now we are ready to use the previous analysis for the evaluation of the hole spectral 
function. In fact, the spinon contribution to the Green function Z p (Q,t) appearing in ([[]) 
can be calculated within the described formalism leading to the expression: 

Z p {R,t) =< q Q \e iipR - Hjt+Eot) \y > (11) 

The dependence on the total momentum p occurs only through the phase shifts (0) defining 
the boundary condition to Hj. The asymptotic behavior of Z p (R,t) can be analytically 
evaluated as: 

Z P (R, t) ~ ^ _ Vst)Xv { P )+^ R + Vat )x v (p)-A ( 12 ) 
showing that singularities characterized by different exponents (|9|) occur at wavevectors 
Q v = (2v + 1)Qf- The additional critical exponent A can be expressed in terms of the 
phase shifts: A = (5r + 5l + 27rz/)(5 R — 5l)/{2it) 2 . When the Fourier transform of Eq. ([12]) 
is substituted into the asymptotic form of the Green function (0) we obtain the anticipated 
expression (|3]) which constitutes the main result of this Letter together with the analytical 
evaluation of the critical exponent x v (p) in Bethe ansatz soluble models. Due to the allowed 
values for v in (|) singularities at Q v = (2v + 1)Qf are predicted, determining "shadow 
bands" in the spectral function. 

From a physical point of view, our results show that the spectral function of Mott 
insulators and superconductors is characterized by branch cut singularities with exponents 
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depending, in the latter case, on the momentum p of the injected particle. This feature 
is shared by all the models we have investigated: The Bethe ansatz solvable Hubbard and 
tJ models and the J — > limit of the Uxy model [[HJ. We believe that it is a general 
feature of hole motion in ID correlated systems thereby providing definite predictions for 
the analysis of photoemission experiments in quasi one dimensional systems characterized 
by a gap either in the charge or in the spin spectrum. 

It is a pleasure to thank M. Fabrizio for extensive discussions and to acknowledge kind 
hospitality at SISSA (AP), Cantoblanco University (SS) and ISI foundation (EU Contract 
ERBCHRX-CT920020). 
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FIGURES 

FIG. 1. Correlation exponent (or dressed charge) K p for the Heisenberg model as a function of 
magnetization. The line is obtained by numerical solution of the integral equation for £22 of Ref. 
§• 

FIG. 2. Lowest critical exponent X u (p) as a function of the total momentum p of the spin down 
hole in the t J model at J = 2 with several magnetizations m. Results are obtained by numerical 
solution of the set of integral equations defining the correction to scaling of the one hole ground 
state energy jy| 

FIG. 3. Comparison between the analytical results of Fig. 2 for m = ±0.25 and Lanczos 
diagonalization. In the latter case, the exponent X v (p) is obtained by a size scaling of the numerical 
evaluation of the overlap ( p as defined by Eq. (|8|). Solid line: analytical results, open dots: Lanczos 
data for m = 0.25, full dots: Lanczos data for m = —0.25. 
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